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Abstract 

oo " 

-^ ' Tensor products of Af random unitary matrices of size N from the circular unitary 

ensemble are investigated. We show that the spectral statistics of the tensor product of 
random matrices becomes Poissonian if M = 2, iV become large or M become large and 
N = 2. 
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H ■ 1 Introduction 



Random matrices proved their usefulness in describing the spectra of quantum systems, the 
classical analogues of which are chaotic O [11]. In particular, spectral properties of the evo- 
lution operator of a deterministic quantum chaotic system seem to coincide with predictions 
obtained for the circular ensembles of random unitary matrices. The symmetry properties 
of the system determine which ensemble of random matrices is applicable. Specifically, if 
the physical system does not possess any time-reversal symmetry, one uses random unitary 
matrices of the circular unitary ensemble (CUE). 

Statistical predictions obtained for ensembles of random matrices are also useful in ana- 
lyzing generic properties of entangled states |12[ [7] . In the theory of quantum information 
one deals with composite quantum system described in a Hilbert space with a tensor product 



structure. Thus any local unitary dynamics can be represented as a tensor product of unitary 
matrices, each describing time evolution of an individual subsystem. 

Consider a quantum system consisting of M subsystems. For simplicity we shall assume 
that each of them is described in N dimensional Hilbcrt space, so that any local unitary 
dynamics can be written as, U = Ui ® ■ ■ ■ ® Um, where Uj € U{N) ioi j = 1, . . . , M . If the 
unitary dynamics of each subsystem is generic, the matrices Uj can be represented by random 
matrices from the CUE. 

The main aim of the present work is to analyze properties of the tensor product of random 
unitary matrices. We show that when either A^ = 2 in the limit of a large number M 
of subsystems, or when Af = 2 in the limit of large subsystem size N , the point process 
obtained from the spectrum of C/, properly rescaled, becomes Poissonian, in the sense that its 
correlation functions converge to that of a Poisson process. 

This paper is organized as follows. In section [2] we provide some definitions and introduce 
our main results, Theorem[T]and[21 and their corollaries; we also provide numerical simulations 
that confirm the results. Section [3] provides the proof of Theorem [T] and of Corollary [TJ while 
Section m is devoted to the proof of Theorem [5] and of Corollary [21 

2 Spectral statistics for tensor products of random uni- 
tary matrices 

The spectral statistics for two ensembles of unitary matrices will be the focal points of our 
investigation. The first case involves two unitary N x N matrices, whereas in the second 
we consider the tensor product of M two-dimensional unitary matrices. As usual, we are 
interested in spectral properties in the asymptotic limits of large matrices, i.e., respectively, 

iV — >• oo and M ^ oo. 

2.1 Background and basic definitions 

We recall some standard definitions and properties of some ensembles of random unitary 
matrices. The simplest situation is a diagonal unitary matrix with eigenvalues being indepen- 
dently drawn points on the unit circle. Such matrices form the circular Poisson ensemble, 
CPE for short. The name reflects the fact that for large matrices the number n of eigenvalues 
inside an interval of the length L << 27r is approximately Poisson-distributcd 

P(^'") -I 

with parameter A = N/2tt. 

Our main interest will be in unitary matrices of size N x N drawn according to the Haar 
measure on the unitary group U{N); such a matrix is called a matrix from the CUEn, where 
CUE stands for circular unitary ensemble. 

Let ^Af be a CUEjv matrix. Denote by (e'^j )jLi its eigenvalues, where we assume that 
the eigenphases 9^ belong to the interval [0, 27r). The random vector {9^ , . . . , 9^) possesses 
a density PctjEn with respect to the Lebesgue measure, which was given by Dyson in his 



seminal paper [3], 

PcuE„«,...,<) = C^ n Ie^'"-e<r (2.1) 

l<fe</<Af 

This expression can be rewritten in the fohowing form (sec Paragraph 11.1 in [9]) 

,7V 



PcuE„«,...,O=Cw(27r)^det[5w(0f-0r)];\^, 



where 



1 sin — 

Sn{x) = ^. (2.2) 

2tt sm ^ 



In particular 

TV 

^^(0) = ^- 

Zn 

The set of cigcnphases of a random unitary matrix can be seen as an example of a point 
process xn on the interval [0, 27r) related to these eigenphases, by which we mean a random 
collection of points {9^ , . . . , 9^} or, in other words, an integer-valued random measure 

N 
k=l 

where Ix denotes the indicator function of X . 

A possible way to describe a point process is to give its so-called joint intensities or, as 
physicists usually say, correlation functions p^ : (IR+)'^ — > K.+ , fc = 1,2, . . .. In our case 
they might be defined simply as 

p^(xi, . . . ,Xfc) = lim -^P (3ji, . . . ,j, \9f^-xs\<s, s^l,...,k). (2.3) 

It is known [3] that the process xn is determinantal with joint intensities 

Pk{xi,...,Xk) = det[SN{xs -a::t)]s,t=i ■ (2-4) 

(Recall that a point process is called determinantal with kernel K if its joint intensities 
can be written as pk{xi, . . . ,Xk) = det[K{xi,Xj)]'l;j^i.) For CUEjv matrices, due to the 
translation invariancc of the measures we have that K^ixi, Xj) = K^ixi — Xj), hence a kernel 
is given by a function Kn{x) of a single variable. We refer to [1] for more background on such 
determinantal processes. 

By definition, the joint intensity p^ equals N\/{N ~ k)l times the k dimensional marginal 



distribution of the vector {9^ , . . . ,9^). Thus 



^^ ^'^' ' det[5^(x,-x,)]:,^i = l. (2.5) 



If we rescale properly the eigenphases of a CUEat matrix it turns out that they exhibit 
nice asymptotic behavior. Namely, it is clear that the point process {-^9^ , . . . , ^9^} is 



determinantal with the kernel ^Sn ( j^%^ ) • Thanks to the fact that this function converges 
when A^ — > oo, we can give a precise analytic description of the limit of the probability 
P {£ef i A,..., ^Q^ i A) , where A C K+ is a compact set (see Theorem 3.1.1 in J]). 

In the case of CPE matrices the situation is even simpler. The point process related to the 
rescaled (by the factor j-) eigcnphases of a CPEat matrix behaves for large A'^ as a Poisson 
point process with the parameter A = 1. 

For point processes, related to the correlation functions is the notion of level spacing 
distribution, denoted by P{s), which is defined for a point process {a^i, . . . ja^N} of the 
properly rescaled eigcnphases {'&j)f^i of a random A^-dimensional unitary matrix by 

1 1 ^ 

where 

si=a {^[ +2tt- d'j^) , sj = a{^'^ - 79;„i), K .7 < N, (2.7) 

and {'&'j)f^i is the non-decreasing rearrangement of the sequence {'&j)f^i. The scaling factor 
a is chosen so that the mean distance Es-,- between two consecutive rescaled eigcnphases is 1. 
In the cases of a CUEat or CPE^r matrix, one has ct = ■^. We should bear in mind that the 
level spacing distribution of the Poisson point process with the parameter A = 1 is exponential 
with the density 

P{s)^e-^. (2.8) 



Moreover, it is easy to check that 



-PcPE„ (s) — > e 



Of course, the limit for the CUEtv is different. 

2.2 Statement of results 

Wc now present our main results for the two cases under consideration. 

2.2.1 M = 2, TV large 

We begin by considering two independent CUE matrices A and B of size N. We are interested 
in the asymptotic behavior of the eigcnphases of the tensor product A(E} B. Our first main 
result is the following. 

Theorem 1. Let {9j)jLi and {(j>j)^^i be the eigcnphases of two independent CUE^r matrices 
A and B. Define the point process gm of rescaled eigcnphases of the matrix A® B as 

N 

aN{D) := Y, l{Nl(e,+^, mod 2.)eD)' ^"'^ ""^^ ^''"^^"^^ '^^ ^ ^ ^+- (^.9) 

k.l = l *■ '" -• 

Let p^ , fc = 1, 2, . . . be the intensities of the process a^ ■ Then 

Pk > 1, (2.10) 



uniformly on any compact subset o/(M+) 



k 



Thus, Theorem [T] relates the statistical properties of a properly rescaled phase-spectrum 
of a large CUE at (g) CUEat matrix to those of a Poisson point process. A (not immediate) 
corollary of the convergence of intensities is the following. 



(2.11) 



Corollary 1. For the point process a^ defined in (|2.9[) . 

P {aN has no rescaled cigcnphase in the interval [0, s]) 



'(aw([0,.s]) = 0) >e- 



Af-i-oo 



s > 0. 



In particular 



^cuEjv «i cuEjv (s) > e 



TV-i-oo 



(2.12) 



where the level spacing distribution PcuEn «> cuEn (*) ^s defined by 




Figure 2.1: The level spacing distributions P(s) for the tensor products of random unitary 
matrices CUEn «) CUEat for iV = 2 (v), iV = 3 (D), N = 20 (o). The symbols denote the 
numerical results respectively obtained for 2^^, 2^^, 2^'^ independent matrices, while the solid 
line represents the exponential distribution (|2.8p . 



Our numerical results support ()2.10|) . i.e. the level spacing distribution of the tensor 
product of two random unitary matrices of size N is described asymptotically by the Poisson 
ensemble. The numerical data presented in Figure 12.11 reveals that just for iV = 20 the 
difference between Pcue„ «> cue„ («) and Pcpe„ («) is negligible. 



2.2.2 iV = 2, M large 

We next consider M independent CUE2 matrices Ai , . . . , Am and study the asymptotic prop- 
erties of the phase-spectrum of a matrix Ai ^ . . . (E) Am- Our main result is as follows. 



Theorem 2. Let 9^9'^, j = 1,...,M be the eigenphases of independent CUE2 matrices 
Ai, . . . , Am ■ Define the point process tm of the reseated eigenphases of a matrix Ai<^ . . .® Am 
as 

tm{D):= Y. ^{--^{el^+...+9lY mod2.)eDy for any compact set D (l^+. 

e=(£i,...,eM)e{l,2}A^ ^ " J 

(2.13) 
Then, for each k there exists a continuous function 5k '■ M+ — ^ K+ with dk(P) = so that for 
any mutually disjoint intervals Ii, . . . , Ik C M-f 

r ¥{TM{Ii)>0,...,TM{h)>0) ^,, , , , ,,,.. 
hmsup — — < (1 + 4 (max \Ij\)) , 

M->oo \tl\ ■ ■ ■ ■ ■ \-lk\ 3 



'(rM(/i)>0,...,TM(/fc)>0) 

"m^'s; iJii • ...• i/fci J 



liminf .^ , '"'1^1 ' > (1 - 4(max |/j|)) 



Note that the statement of Theorem[2]is weaker than that of Thcorcm[T] This is due to the 
fact that stronger correlations exist in the point process tj\/ , which prevent us from discussing 
the convergence of its intensities to those of a Poisson process. The mode of convergence is 
however strong enough to deduce interesting information, including the weak convergence of 
the processes. We exhibit this by considering the behavior of the level spacings when M tends 
to infinity. 

Corollary 2. For the point process ti\i defined in (|2.13p we have 

P (tm has no eigenphase in the interval [0, s]) 

= P(rM([0,,s]) = 0)-^ ^e-^ s>0. ^^'^^^ 



M- 



^cuEr(s)T^ ^e-^ (2.15) 



In particular 



where the level spacing distribution P„yj-^isM{s) is defined by (j2.6[) . 

The relevant numerical results which confirm (|2.15p are shown in Figure 12.21 Again we 
may observe that it is enough to take relatively small M in order to get a good approximation 
of the spectrum of a matrix CUEf" by the Poisson ensemble. 

2.3 Discussion 

The convergence exhibited in Theorems [T] and [21 and in their corollaries, is arguably not 
surprising: taking the tensor product introduces so many eigenphases (iV^ in the case of 
Theorem[Tl 2*^ in the case of Theorem [2|) that, after appropriate scaling, the local correlations 
between adjacent eigenphases are not infiuenced by the long range correlation that is present 
due to the tensorization. One should however be careful in carying this heuristic too far: 
well known superposition and interpolation relations, see [3] and the discussion in [U Section 
2.5.5], show that the point process obtained by the union of eigenvalues of, say, a GOE^r and 
GOEjv+i independent matrices, is closely related to that obtained from of a GUEat matrix, 
and thus definitely not Poissonian. This phenomenon had been also discussed in the physics 




Figure 2.2: Level spacing distributions P{s) for the tensor products of random unitary matri- 
ces CUEf^' for M = 2 (V), M = 3 (D), M = 8 (o). The symbols denote the numerical results 
respectively obtained for 2^'^ , 2^^, 2^^ independent matrices, while the solid line represents the 
exponential distribution p.Sp . 



literature |10| . Compared to that, the tensorization operation appears to strongly dccorrclate 
eigenphases on the local level. 

It is natural to try to generalize Theorems [1] and [2] to other situations, where either N or 
M are finite but not necessarily equal to 2, or both N and M go to infinity. While we expect 
similar methods to apply and yield similar decorrelation results, there are several technical 
issues to control, and wc do not discuss such extensions here. 



3 Tensor product of two N x N unitary matrices 

We prove in this section Theorem [T] and Corollary (TJ that correspond to the case M — 2 and 
N large. We start with an elementary observation. Recall the kernel Sn, see 



Lemma 1. For any N > 1 



N 
sup|S'Ar(a::)| = — . 



(3.1) 



Proof. Wc show inductively that 

I sin(nu)| < n\ smu\, for m g R,n > 1. 

Hence 



\Sn{x)\ = ^ 



sm 



(^f) 



< 



N 
2^' 



For a; = we have equality, which completes the proof. D 

Proof of Theorem]^ We begin with setting Si, . . . , Xfe > and recalling that by definition 



7V2 



7^3 (6*2;^ + 0j^ mod 27r) € (i^ — £, is + e) 



Let us first of all get rid of the addition modulo 27r noticing that the event, probability of 
which we want to compute, is the sum of 2'^ mutually exclusive events occurring when Oi^ +4'jb 
is in the interval [0, 2tt) or [2tt, Att). Thus we can write the sought after probability as 

^ v(3^'^ys0,^+(j),^eiT],-2n + x,-e,Ji,-2TT + x,+e)Y (3.2) 

where we denote Xs = -^is and e = ;^e- Let us now concentrate solely on the first term 
corresponding to the index (ry) = (ryi, . . . , ry^) = (0, . . . , 0) (the other terms can be dealt with 
in the same manner). In order to take an advantage of the independence we observe that the 
considered quantity equals 

if^oo ^-^ _ \ (j) 0j^ e {xs - 27r4/X - e, Xs - 27r4/i^ + e)) ' 

1tvIb)k<Xb 

where the constrains 2'Kts/ K < Xg are the result of the fact that 6'j^ + (pj^ e (0 • 27r + Xg — e, • 
2Tr + Xs + s), for (77) = 0, so, in particular, that 9i^ < Xs + e. Exploiting the independence we 
obtain that the last expression equals 

K 

Y^ P(3(i)V.s 0i^ e{2TTis/K ~n/K,2TTl,/K + n/K)) 
^i,...A.=i (3.3) 

2Trf„/K<XB 

■ P i3ij)ys 0J-, e {xs - 2tt£,/K -e,Xs- 2-k(.s/K + e)) . 

Now observe that for a dcterminantal point process {ctj}f^i with a kernel K and fixed numbers 
Ml, . . . , Ufe we have 

P (3(i) G {1, . . . , A^}*^ V.S == 1, . . . , fc a^^ e (ws - (5, u, + S)) 

= E E A,(Mi,...,Mfc)((2<5)Pdet[if(u,(s,i),M,(M))]s,,=i+o(<5^)). ^^'^^ 

p=i 7ree(fc,p) 

where G{k,p) is the collection of all partitions into p non-empty pairwise disjoint subsets of 
the set {1, . . . , fc}. By this we mean that if n is such a partition then 

TT = {{7r(l, 1), . . . , ^(1, tt7r(l))}, . . . , {^(p, 1), . . . , 7r(p, tt^(p))}}. 



where tt7r(g) is cardinality of the q-th block of the partition tt. Moreover, to conipactify the 
notation, we attach to a partition tt the function A^r : M'"' — >■ {0, 1}, defined as 

AttIUi, . . . , Uk) ~ 1{m^(i i)=...=«^(i ,^(i)),...,M^(p_i)=...=u^(j, ,^(p))}(Ul, ■ • ■ , Uk). 



Applying this to formula 

K k 



we obtain 



E E E A., ((27r^./i^)ti) A., ((x, - 27r^,/i^)ti) 

ii,...,ik=i pi,p2=iiTiee{k,pi) 

2Tri,/K<x, TT2e&{k,p2) 

(2_\ Pi 
— j det [SN{27r£^,is,i)/K-27T£^,^t,i)/K)Y;,^^+o{l/Kn 

■ ((2er det [Sn K,(s,i) " ^^^Ms-d/K " ^-^(ta) + 2^4,(*,i)/if)]',^i + o{eP')) . 

Performing the limit i^ — > oo we notice that only the terms corresponding to p2 ~ k do not 
vanish, for, otherwise, A^j would give nontrivial relations for {£) which altogether with A^ri 
make the sum over {£) of at most 0(Xpi~^) terms. Recall that e/e = 2tt/N'^ . Thus, taking 
the limit e — )- 0, the extra factor {2'k/N'^Y is produced, so we finally find that the considered 
term contributes 



E]^ E ji;^pjy,a^)My^^---^v^)'^^' 



-jrSNiyTrisS) - VTritS)) 



P 

s,t=l 



•det 



— Sn{xs ~ys-xt +yt) 



d'Hpiyi,...,yk) 



s,t=l 



to p-^{xi, . . . , Xk), where Up denotes the p-dimensional Hausdorff measure in R*^. As already 
mentioned the other terms in (j3.2p can be calculated in a similar way, only the limits of the 
integration have to be changed. Summing up, we get 



N / 

Pk yxi, 



,Xk) 



Z^ 1^ Nf^-P ^ {2n)P Ja 



('7)6{0,1}''P=1 



TTee{k,p) ^ ""'^ •^'^("i 



X^{yi,...,yk) 



•det 



•det 



^'5'Ar(y7r(s.l) - J/7r(t,l)) 



s,t=l 



27r 



SNi2mi, + Xs -ys - ^Trrjt -Xt + yt) 



|dHp(yi 



,yk), 

(3.5) 



where the subset yl(j,) of [0, 2tt)'' is the set of all (j/i, . . . , j/fe) such that cither ys < Xg iiris = 0, 
or j/s 1^ Xs ii rjg — 1 for s = 1, . . . , fc. 

To proceed we have to investigate the asymptotic behavior of the integrand in p.Sp . We 
will do it again only for (77) ~ (0, ...,0), observing that an adaptation to other terms is 
straightforward. We start with the term p = k. Then the integrand is a product of two 



determinants of matrices of size k, so applying to each of them the permutation definition 
and extracting the term referring to the trivial permutations, we find it equals 



N 



5jv(0)j + ^ Sgnaiig'a.T\\_—SN{yi-ya{i))W_^SN{Xj-yj-Xr{j)+yr{3)), 

(3.6) 
where the summation involves all permutations a and r of k indices. The first term 
{2ttSn{0)/N)'^'' = 1, after substituting in p.Sp . gives simply 



(';)6{o,i}'= <"' ^ 



2k 



We will show that the second term in p.6p after being put into ()3.5p vanishes in the limit. 
We consider here only the case A: = 2 to explain the main idea. The terms involving more 
factors can be treated along the same lines. The sum over a and t reduces to 



^'5'Ar(0) j I ( -j;^SN{yi - y2. 



—Sn{xi -yi- X2 + 2/2) 



SNiyi -2/2) 



(3.7) 



SNixi ~yi - 2:2 + ^2) 



Let us for instance deal with the last term in equation p. 71) . Putting it into p.5p we arrive 
at 



i^^LS^'''^''-''^) (^ 



SNixi -2/1 - X2 +2/2) 



Taking a quick look at the integrand we see that the above expression goes to when N ^ 00 

by Lcbesgue's dominated convergence theorem, for j^Sj\r{u) > 0, when u ^ 0, and the 

appropriate bound p.ip follows from Lemma [1] 

For the terms corresponding to fc < p, we easily notice that thanks to the factor „fc_p 
they converge to 0. The proof is now complete. D 

Remark 1. By virtue of formula p.Sp the joint intensities p^ can be estimated as 



sup|pf|<— ^ sup det 

R*" -'^ iii,...,MfceR 



2tt 



SNiUs -ut) 



y'tte(fc,p) / det[S'w(2/s-2/t)]st=id2/i---d2/p, 

p=l J[G,2-n)P 



where '^X denotes cardinality of a set X. Using Hadamard's inequality (see, e.g. (3.4.6) in 
[T]) for the first term, the observation (|2.5p for the second one, and finally p.ip we obtain 



^N\ 



1 



sup|p,|<^(sup 



27r 



^'''E«6(fc'^)n^-^'''^E«6(^'^) 



iV! 



p=i 



(A^-fc)! 



p=i 



[N -k)V 



10 



Due to the well-known combinatorial fact that 

k 

Y^ ^e{k,p)x{x - 1) • . . . (a; - p + 1) = x^ 
p=i 

(tt(5(fc,p) is the Stirling number of the second kind, consult e.g. [5]) we may conclude with an 
useful bound 

suplpfl <A:^/^ □ (3.8) 

Proof of Corollary [J^ For the proof of (|2.1ip we have to calculate the probability of the event 
that there is no rescaled cigenphase in a given interval. This is done in the following lemma. 

Lemma 2. Let x be a point process related to the eigenphases, possibly rescaled, of a CUEjsi 
matrix Am with the joint intensities pk, /c = 1, 2, . . . (so pe = 0, for £ > N ). Then for any 
compact set D 

F{x{D) = 0) = l + J2^—L p,. (3.9) 

t=i ' •^■°* 

Proof. Clearly, we have 

TV 



'(x(z?)=o) = i-^p(xp) = fc). 



k=l 



One way to compute the probability P {x{D) = k) is to use the notion of Janossy densities 
JD,k{xi, ■ ■ ■ ,Xk) (see Definition 4.2.7 in [1]). They can be expressed in terms of the joint 
intensities as 

JD,k{xi,...,Xk) = Y —{-ly pk+r{xi, . ■ . , Xk, D, . . . , D), (3.10) 

~ r 

where 

Pk+r{xi,...,Xk,D,...,D) = / pk+r{xi,...,Xk,yi,...,yr)dyi-- -dyr. (3.11) 

r 

They exist whenever 

yf fc:Pfc(^2j_lllMd^,...d;^, <oo, (3.12) 

k-^D" 

which is clearly fulfilled in our case, as p^ = for £ > iV. Moreover, the vanishing of pe 
for large enough i makes every sum in the following finite so we will not have troubles with 
interchanging the order of summations. 

In terms of the Janossy intensities, the probability ¥ {x{D) = k) reads as (see Equation 
(4.2.7) of m 

F{x{D) = k) = ^ [ jD.k{xi,---,Xk)dxi---dxk, (3.13) 

fc! Jjjk 



11 



and, consequently, 

N 

P(X{D) = 0) = 1-J2-J^ 



iD.k 



k=l 



^ fc! .Id" ^ H 



fe=i 



/0/c+r(a;i , . . . , Xfc, Z?, . . . , i:))d.Ti . . . dxk 



1-EE 



fc>l r>0 



>0 



Pk+ 



r=^-j:j: 



1 (-1) 



i-k 



k>ie>k 



k\ {(.-k)\ Jai 



1-E 



£>1 



Lfe>l 



E M 






P£ 



-1 + E 



f>i 



Ml! 



Pf- 



D 



Lemma [2] applied to the process a^ yields 

P(a^([0,,s])=0) = l + ^^/ 



i>i 



N 
Pi ■ 



To pass to the limit A^ — >■ oo we need an appropriate bound on the intensities p^ . In Remark 
[T]we showed that \p^\ < tl"^ (see p.Sp ). Therefore, by Lebesgue's dominated convergence 
theorem, we get 



lim P {a^{% s]) = 0) = 1 + 5] ^^ / lim pf = 1 + ^ 






This completes the proof of (|2.1ip . 

The formula (|2.12p follows now from a relation connecting the probability E(0, s) that a 
randomly chosen interval of length s is free from eigenphases with the level spacing distribution 
P{s), daSl) (see equation (6.1.16a) in [9]), 



We have just showed that hmAr^oo P (crjv([0, s]) = 0) = E{0; s) = e^'^. Thus, indeed 



(3.14) 



lim PcuEjv «) CUE„ (s) = --^e '^ = e "*. 

Af-i-oo ds 



n 



4 Tensor product of M unitary matrices of size 2x2 

Now we will prove Theorem [2j In the course of the proof we will need three lemmas. Let us 
start with them. 



12 



Lemma 3. Fix a positive integer s and a number 7 S (0, 1/s). For each positive integer n 
let us define the set £„ = {i = {ii, ■ ■ ■ ,£s) \ "^ 3 £j > 0, J2'j=i ^j = "-}• Then 



El n! Y^ 1 ^^• 

Here, we adopt the convention that (.1 ^ iil ■ . . . ■ £s^-- 
Proof. First observe that 



^0. 



(4.1) 



[fn] 



El n! Y^ 1 "■■ Y^^ 1 



e,3j i,/n<-r 



£,£i/n<7 

L7"J 






1=0 



n 

s" V n - £1 



k—7i— [7nJ 



E 



(»-4)! 



1-i 

s 



Let Xi,X2,..- be i.i.d. Bernoulh random variables such that P (Xi = 0) = l/s = 1 — 
P {Xi = 1). Denote 5'„ = Xi + . . . + X„. Then the last expression equals sP (5„ > n — [7?iJ ) 
and can be estimated from above as follows 



sP {Sn>n- 771) 



'^" ^^" > i - 7 ) < sexp ( - 2n(l/s - 7)^^ 



^0, 



where the last inequality follows for instance from Hoeffding's inequality (see [8]). 



n 



Lemma 4. Let X be a random vector in M" with a bounded density. Let A : M" — > M.'^ 
be a linear mapping of rank r. Then there exists a constant C such that for any intervals 
/i , . . . , /fe C M 0/ finite length we have 

P(AXG/iX...4)<C|/,J-...-|/,J, 

where I < ii < . . . < ir < k are indices of those rows of the matrix A which are linearly 
independent. 

Proof. Let ai,. . . ,ak G K" be rows of the matrix A. We know there are r of them, say 
ai,...,ar; which arc linearly independent. Thus there exists an invertible r x r matrix U 
such that 



E. 



where e^ € M" is the z-th vector of the standard basis of R" . Notice that 

V{AX e /i X ... x/fe) <V{U-^EX eLi X ... x/,,) =V{{Xi,...,Xr) e C/(/i X ...Lr)) 
< C|C/(/i X . . . X /fe)| = C| det C/| • |/i| • . . . • |/,|, 

for the vector (Xi, . . . , Xr) also has a bounded density on W . This finishes the proof. D 





'ai' 




'ei' 


u 


'. 


= 






Or 




.'^■r. 
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Lemma 5. Let A be a matrix of dimension k x j , with entries in {0, 1}, and satisfying the 
following conditions 

(i) no two columns are equal. 

(ii) no two rows are equal. 

(Hi) one column consists of all Is. 

Then, the rank of A is at least niin(fc, [logj j J + 1). 

Proof. (Due to Dima Gourevitch) Denote r — rankyl. The assertion of the lemma is equivalent 
to the statement that 2'' > j and if 2^ — j then r = k. 

We may assume without loss of generality that the first r rows of A are linearly independent 
and the others are their linear combinations. Under this assumption, if two columns are 
identical in the first r coordinates then they are identical in all coordinates. By condition (i), 
such columns do not exist. Therefore the r x j submatrix B which consists of the first r rows 
has distinct columns. As a result j < T' . 

Now suppose i — T" . If fc > r, consider the r + 1 row of A. It is a linear combination 
of the first r rows. Since the columns of B include the column e^ = (0, .., 0, 1, 0, .., 0) for all 
i = 1, . . . , r, the coefficient of each row is either or 1. Since A includes a column of Is, the 
coefficient of exactly one row is 1, and all other coefficients vanish. Thus, the r + 1-th row is 
identical to one of the first r rows - in contradiction to condition (ii). D 

Proof of Theorem [H Fix an integer k > 1 and finite intervals /i , . . . , /^ C K+ which are 
mutually disjoint. We need to compute the probability of the event {tm{Ij) > 0, j = 1, . . . , fc} 
which means that in each interval Ij there is a rescaled eigenphase. Such eigenphase is of the 
form ^ {91' + ... + e\\' mod 27r) for some e = (ei, . . . , em) € {1, 2}"'. Therefore 



where 



{TM(/j)>0,j = l,...,fc} = |jA„ 



{^ e^ 2-K 1 

^^ = I E^^' '"od 271 e ^/„.7 = 1, . . . , fc j, 



(4.2) 



J, 



and e runs over the set 

f - {[e^&;;:M I ^ ^ {1,2}, e^^e\ioTu^v,u,v^l,...,k] (4.3) 

of all k X M matrices with entries 1, 2 which have pairwise distinct rows e^ = (ej, . . . , e'j^j) e 
{1, 2}^^, j = 1, . . . , fc (j'-th row e^ describes the j'-th eigenphase and since intervals are disjoint 
we assume the rows are distinct). Column vectors are denoted by ei ~ [e^, . . . , e*^]"^, i = 
1,...,M. 

We say that e is had if the collection of its vector columns {ei,i = 1, . . . , M} is less than 
2^ . Otherwise e is called good. 
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Obviously, 




The strategy is to show that the contribution of bad e's vanishes for large M while good e's 
essentially provide the desired result Y\n I -01 when M goes to infinity. So the proof will be 
divided into several parts. 

Good e's. The goal here is to prove 



lim lini 



1 



maxj |/j|->.0 M-i-oo \Ii\ • . . . • |/j. 



u ^^ 



1, 



i good e's 

with the required uniformity in the choice of the disjoint intervals Ij . By virtue of 



(4.4) 



^ P(A,)- Y. ]P(^.nA^)<P \J aA< Y, P(^) 



good e's 



it suffices to prove that 



good e, € 



k good e's 



good e's 



lim > ] 

eood e's 



'(^e)=ni^.- 



(4.5) 



uniformly, and that the correlations between two different good epsilons docs not matter 

1 X- 



lim sup lim sup y-r , ^ , / 

max, |7,|-).0 M^oo \\.\h\ g^^^ ^ ~ 



P {A, n Ai) = 0. 



(4.6) 



Let us now prove (j4.5p . The proof of (|4.6p is deferred to the very end as we will need the 
ideas developed here as well as in the part devoted to bad e's. 

Given e € £ and a vector a — [ai . . . afe]^ £ {1, 2}*^ we count how many column vectors 
of e equals a and call this number £„. Then X]q,-^q = M. Note that e is good iff all €qS 
are nonzero. The crucial observation is that the probability of the event A^ does depend 
only on the vector i = {^a)ae{i,2}'' associated to e as described before. Indeed, the sum 

where 



"^ mod 27r is identically distributed as the random vector J^a '/'(o^i ^a) mod 2tt, 



■4'{a,ia) 



ijjkiaja) 



41 



41 . 



gai 



■ 4«a 



mod 27r 



(4.7) 



is a sum modulo 27r of i.i.d. vectors. Note that the distribution of ip{a,ia) does not depend 
on the choice of indices ii, . . . , i^^ but only on a and £a- Consequently, denoting by £i the 
set of all e's such that there are exactly £„ indices 1 < ii < . . . < i^^ < M for which 
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Eij = . . . = Ei^ = a, we have that the value of P(^e) is the same for all e G £i. Clearly 
tJ££ = -^, whence 

Y, P{Ae)^ Y. TT^ H ^(a,£„)niod2^e J1X...X jJ . (4.8) 

good e's good ts \a£{l,2}'' j 

The idea is to identify those terms which will sum up to Jl 1-^*1 ^'id the rest which will be 
neglected in the limit of large M . To do this, set a positive parameter 7 < 1/2'^' and let us 
call a good i very good {v.g. for short) if ia > iM for every a and quite good [q.g. for short) 
otherwise. We claim that 

^ V(a,^Q) mod 27r e Ji X . . . X Jfcj < C" JJ |J,|, for a good i, (CI) 



and 



(^^(a,£„) mod 2^ G Ji X . . . X Jfe) = H]^ A + -^"j , |r,| < C 



(C2) 
/or a very good £, 

where C is a constant (from now on in this proof we adopt the convention that C is a constant 
depending only on k which may differ from line to line). 

Let us postpone the proofs and sec how to conclude (|4.5p . Notice that UK^ = 2m- H 1^1- 
Thus applying (|Cip we obtain 



q.g. i's q.g. ^'s ■" 

By Lemnia[3]it vanishes when AI — > 00. Now we deal with very good fs writing with the aid 
of (lC2l) that 

Y P{Y^{a,e^) mod2n e J,x ...X J,) =l[\lJ ^ ^^ 

V.g. £'s \ v.g. ts 

v^ 1 Ml re \ 



v.g. £ s / 

The first term in the bracket approaches 1 in the limit M — > 00 due to Lemma [3l while the 
second one approaches as it is bounded above by C^fj. 

Proof of (|C1[) . Let us define the vectors 

eJ = (2__J,l,2__^)e{l,2}^ j = l,...,fc. 

Since £ is good, in particular we have that £e- > 0, so denoting the random vector ip(ej,£e) 
by *J we have EaV'(a,^a) = (*^ + • • • + *^') + Ea^{ei,...,efc} V^("'^")- By independence 'it 
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is enough to show that the random vector v]/ = ij/i + . . . + v]/*^ mod 27r has a bounded density 
on [0, 27r)''. Equation (|47)) yields that 

vi/^ = (r„...,r„x„r„...,r,), 

"^ V ' ^^ V ' 

j-1 k-j 

where {Xj,Yj) are independent random vectors on [0,27r)^ with the same distributions as 
the vectors {6l + . . . + 9j mod 2tt,91 + . . . + 6"^ mod 27r) respectively. Clearly, the vector 

{Xj, Yj) has a bounded density on [0, 27r)^ because the vector {6l,6l) has a bounded density. 
Therefore the vector {Xi, Yi, . . . , Xk, Yk) has a bounded density on [0, 2tt)'^'^. A certain linear 
transformation with determinant 1 maps this vector to (^^ + ... + '^'',Yi, . . . ,Yh) which 
consequently also has a bounded density. One projects it to the first k coordinates and then 
takes care of addition modulo 2tt obtaining that ^ has a bounded density, which finishes the 
proof. n 

Proof of (jC2[) . Given a vector a E {1,2}'^ let 9" denote the random vector in [0,2tt)'^ 
identically distributed as the vector {6"^ , ■ ■ ■ , 0"''). Take its independent copies Qf, 83 , . • . 
such that the family {O", O2 , • . ■}ae{i,2}'' ^l^o consist of independent random vectors. Then 

EO" = [tt, ...,7r]^, and 

PiM = P ( XI ^("' ^") '^^"^ 27r e Ji X . . . X Jfc I = P ( XI X] ^° 1^^°^ 27r e Ji X . . . X Jfc j 

M-l / f„ \ 

= X P X£®"^ (^1+2^*1) ^•■■^("^fe + 2^*^-) 

ii,...,ifc=0 \ a 1 = 1 J 

E E ^^Stt^ ^ 7tj('^i + 2-(H - M/2)) X . . . X -=(J, + 2.fc - M/2)) . 

a 1=1 V V V y 

To ease the notation we introduce new indices 



j = 


/ M M\ r A/ M M 1'' 




^J,M = -/^(^l + ^T^Jl) X • • ■ X -7^(-^k + 2TTJk), 


the vector 





s.^-T.t''''^^' 



a 1=1 ^ 

Now we intend to use the local Central Limit Theorem of [5] . Indeed, due to independence 
such a theorem should hopefully yield that Sm has a normal distribution for large M. To be 
more precise, let us consider the matrix Cov Sm ~ X]q m ^^'^ ®" ^'^'^ its eigenvalues. Since 
for any x € M*^ 

x^{GqySm)x = y -^a;^(Cove")a; < max|l Cove^H^/^ j^^p^ 

a ^ ^ ^ 

C 
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it is clear that the largest eigenvalues are uniformly (i.e. with respect to M) bounded by C, 
which depends only on fc. To provide an uniform bound for the smallest eigenvalues let us 
observe that (recall that e^ is the vector (2, . . . , 2, 1, 2, . . . , 2)) 



k p 
X 



{Gov Sm)x > Y, -^x^{CovQ''')x > 7^'^ ( ^ Gov O''^ ) x > 7 • ^|a;p, 

where the second inequality is because £ is very good. 

It is a matter of direct computation to see the last inequality as for k > 2 wc have 
Ej^i Cove-^^ = ((fc - 2)7rV3 - 2)[1 . . . 1]^[1 . . . 1] + diag(2 + 27rV3, . . . ,2 + 27rV3) and for 
fc = 1 the smir equals 7r^/3. Therefore, with the matrix Bm given by 

it holds that 

^\x\<\Bmx\<C\x\. 

Therefore the assumptions of [H Corollary 19.4] are satisfied (for the family of independent 
random vectors {0", 62 , • ■ ■}aeii,2}'')i so the vector Bi\iSm possesses a density qm and 

sup (1 + |a;|'=+2) (qM{x) ~ cj,{x) - -^Pm{x)cJ>{x)] = OiNr'^'^), 

where 4>{x) = }—k e~l'^l '^ is the density of the standard normal distribution in M'"' and Pm 
is a polynomial of degree fc — 1 whose coefficients depends on the cumulants of the vectors 
BmQ"- We may put it differently, i.e. 

QMix) = Hx) + ^(pm{x)Hx) + , ^''[f,+. 



jiy ^ '^' ' i + \x\^^ 



llMix) 

for some functions Jm uniformly bounded supj\^ sup^gjji- |/A/(a^)| = C < 00. Therefore, de- 
noting Lj^M = BmKj^m, 

= ^ P {Sm e A'j- m) = ^ P {BmSm e BmK.m) ^Y. f 



E/ 'f'+^uT.f hM = aM + -^^ 



(4.9) 



-.b 



Let us firstly deal with the error term 6a/. Denoting 

^ \Ji\---.-\Jk\ 

we are to show that 



M- 



\bM\ < Ck. (4.10) 
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To do this we estimate the integrated function 



\hM{x)\ < \PM{xMix) + ^^1^^^, =: hix). 



Then |6m| I^TIij Jl ^' Introduce full boxes 

Fj,M = Bm ( ^([0, 2tt) + 27rji) x . . . x ^([0, 27r) + 27rjk) 



l=([0,2^) + 27rji)x...x^^ 
and observe that 

h = ttt^I-F'j.a/Itt r / h< k\Fj.m\ sup h < k|Fj,a/| sup h. 

Since diamF,j\/ < C '^^X-^ > 0, the sets K^j\/ are pairwise disjoint and sum up to 

Sa/[— TrViV/, TTvM)''", we can infer that the sum ^- \Fj^j^i\snpp. ^^ h converges to J^^ h ~ 
C < oo. Hence, this sum is bounded by C and we get (14.101) . 

Now we handle the main term qm- We prove it equals k up to another error K-fcy. Let 

Aj^M '■ R'^ — > R*^ be the linear isomorphism mapping Fj^m onto Lj^m- It equals Bm-Aj^mB^{^, 
where Aj,m is the hnear mapping transforming the box B^j Fj,M onto the box B^^ Lj^m, 
whence I dctylj_A/| = k. Thus, changing the variable we obtain 

/ (l){x)dx ^ n <PiAj,Mx)dx. 

Notice that Aj^mx is close to x, whenever x S Fj_m, for 

I^j,n2; — x\ < diami^j.M, x e Fj,m- 

Consequently, on -Fj,a/, 4'{Aj^mx) is close to 4>{x). Strictly, we use the mean value theorem 
and get 



4){x)(lx = K / (t){x)<lx + K '^4>v{Vx) ■ {Aj^MX — x)dx, 

Lj.M "^ Fj.M ^ Fj^M 

for some mean points rjx G [x, Aj^mx]. This results in 

aM = 2.^ / (j){x)dx ~ K y^ / (j) + K V(j){rjx) ■ {Aj,mx — x)dx 

j J Lj.M j Jfj,m J\J,f,_m 

= k(1- ^ + y2 \7(l)v{Vx) ■ {Aj^MX- x)dx 

\ Jr''\Bm[—kVM,ttVII)'' , JFj,m 



We are almost done. Clearly cm converges to faster that l/VA/, so \cm\ < C/vM. For 
(Im we use the Schwarz inequality and integrability of |V0(772:)| 

MmI < y] / \'^4'{Vx)\\A-j^MX - x\dx < diamFj^M / |V0(7?3,)|da; < -y=. 

j Jf,,m J\JFj,m VM 



This completes the proof of (jC2p . D 
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We have proved claims (|C1[) and (jC2p . so the proof of the part concernmg good e's is now 
complete. Let us proceed to tackle bad e's. 

Bad e's. The goal here is to show that 

.iTc^fU ^e)=0, (4.11) 

Vbad e's / 

again, with the required uniformity. Obviously it suffices to show that X^bad e'- ^ (^e) ^ 

0. Let Tj be the set of those bad e's for which the cardinality of the set {e^, i = 1, . . . , A/} 
equals j. Observe that %Tj < j^^ . With the aid of Lemma [S] we will show that 

Ve e J^i P{A,) < C ■ 2-A^(i+Li°g2jJ) . 0((max|/,|)i+Li°S2jJ)^ ^hen max|/,-| — > 0. (4.12) 

j j 

This will finish the proof, for 

2''-! 



J2 P (Ae) < C ■ 0(max |/,|) Y, /' • 2-*^(i+Liog, jj) 



bad e's J — 1 



(4.13) 



C ■ 0(max \Ij\) y 2"*^(i+Li°g2jJ-iog2 J) ^ q 






A/- 



For the proof of (|4.12p fix e e J'j . We have seen that 

P {A,) ^ P (y2 ^'("> ^") mod 27r e Ji X . . . X Jk) 

and we know that there are exactly j numbers la which are nonzero, say those which corre- 
spond to vectors a^,...,a^ G {1,2}''^. Denote ^^ = -(/'(a^^ctO; * = 1j---:J ^-^d consider the 
random vector Sj = '^^ + . . . + "^^ in R*'. As in the proof of Claim (JC1[) we observe that 
S'j is a linear image of the vector {Xi, Yi, . . . , Xj, Yj). This mapping is given by the matrix 
A ~ [ciuv] where 

Jl, a»„ = l 
a2,-M-j^^ < = 2' "^^' 

By Lemma |4]wc obtain 

P (S'j mod 27r e Jl X . . . X Jfc) < C max (| J,J • . . . • | Jj^ |) = C • 0(max |/j |) • 2"^^'', (4.14) 




where r = rankA. The number r does not change if we replace the 2i-th column of A with the 
vector e with 1 at each its entry, as the sum of 2z — 1-th and 2i-th columns is just e. Taking 
only the columns 1, 2, 3, 5, ... , 2j — 1 we get the matrix B which has the same rank as A. It 
has j + 1 columns and fulfils the assumptions of Lemma[5l Thus r > min(l + [log2(l + j)J , k) 
and when j < 2^ — \ this minimum equals 1 H- [log2(l + j)\ > 1 + [log2 jj. If j = 2'^ — 1 
in the matrix A there must be two identical columns, one with even, say 2m, and one with 
odd, say 2v — 1 index, which means that the u-th and the w-th column of B add up to e, 
so the w-th column may be erased and the rank of B does not change. Therefore we apply 
the lemma to the matrix B with erased the w-th column which is of size fc x j and get again 
r > nun(l + Llog2 j\^k) ^ 1 + [logj j\ . This completes the proof of (|4.12p . 
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Pairs of good e's, i.e. the proof of (|4.6p . We denote by 8i(e) the random vector 
(6*,' , . . . , 6'j' ). By the definition of A^ we may write 



AnAr^\j2 



■"6,(6)1 , „ ^ Ji X ... X Jfc 1 , . 

^ ;~; mod ZTT e , t > • (4.15 

e,(?) Ji X . . . X Jfc I ^ 



Since the intervals Ju and Jy are disjoint for it 7^ w, we may restrict ourselves to those e and e^ 
for which e" j^e" whenever u y^ v, u,v = I, . . . ,k a.s otherwise the event A^nA-^ is impossible. 
However it might happen that e" = e^. Let us count for how many u's it takes place, i.e. 
given s G {1, . . . , fc} let Ps be the set of all considered unordered pairs {e,e} for which there 
are exactly k — s indices 1 < ui < . . . < Wfc_s < k such that e"^ = e^^ , j ~ 1, . . . ,k ~ s. The 
value s = is excluded as e 7^ ?. We have 

k 

Thus we fix s and prove that limsup„^^^^ |^^.|^o limsup^^^^ 1-|^ E{£,?}gp, ^ (^^ H ^j-) = 0. 
There are two cases. A pair {e, e} G T'is can be good which means tt{[?-], i = ^, ■ ■ ■ ,^l} > 
2k+s^ or, otherwise we call it bad. We obtain a decomposition Vs = V§°°'^ U V^^'^. Now for 
a good pair, applying the reasoning already used for bad e's, i.e. combining lemmas H] and [Sj 
we get the estimate 

V{A, n A,) < qjil • . . . • kfeK^max^ 1^,1)' = ^(I^ (R l^^'l) ("^f^ I^^D'' 
But UTT"^ < ]),Vs < (t) • 2('=+'')^^ so 

limsup lim sup V" P(Aen^j)=:0. 

max,-|/,|-).0 M^oo 1 1 1^' I ^^ j^^^jood 

For a bad pair {e, e} we know that there are k + s different rows and at most 2'^"*"* — 1 
different columns in the matrix [-]. Hence we repeat the argument of the part concerning 
bad e's. Namely, first exactly in the same manner as in that part we use Lemma [5] in order 
to establish an appropriate inequality in the spirit of (|4.12p . Then we follow the estimate of 
(|4.13p and conclude that 

lim V PM, nAj) = 0. 

This finishes the proof of Theorem [21 D 

Proof of Corollary\^ Fix A small so that s/A is an integer and divide the interval [0, s] into 
consecutive intervals of length A, denoted U. Let Zi = TM{Ii) and Zi = l{2.>o}- Of course, 

tm([0, s\) = YlilLi ^i- Our goal is to show that tjv/([0, s\) becomes Poissonian in the limit of 
large M , from which the statement of the corollary follows immediately. 
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The proof of Theorem [2] yields the fohowing facts. There exist a sequence Sm^am with 



Hm sup Hm sup 6M,A,k ~ , 

A-i-O A/-S-00 



and a universal constant C such that the following hold. 

P(Z, ^Z,) < CA\ 



(4.16) 
(4.17) 



\ieJk 



where Jk denotes an arbitrary subset of k distinct integers in {1, . . . , s/A}. 
Indeed, to justify (|4.16p notice that 

P(Z,^Z,) = P (tm(/.) > 2) < P|J A, nA^ , 

where e,e G {1,2}^^ and A^ is the event that there is an eigenphase described by e in the 
interval li (see (|4.2p ). The probability of the event A^ n A-^ can be estimated by C ■ 2~^^^ • 
|/ip = 2~'^^^ ■ CA^. To see this, recall (|4.15p and follow the same argument which led to 
estimate (|4.14p (in this case the relevant matrix has the rank no less than 2). It suffices, as 

fU.^?^£ n Ar) < (2") • 2-2^^ • CA2 < CA2. For (PTTI) . observe that 



E [] Z, = El{2.>o,,ej,} - P(tm(/.0 >0,ie .h) , 



i£Jk 



and apply Theorem [5] (with its uniformity statement). 

Let Y, be i.i.d. Bernoulli random variables with P (Yi == 1) = 1 - P (Yi == 0) = A. By 
(|4.17p we have that for any integer £, 



lim sup lim sup 

A-J-O M-foo 



,s/A \ /s/A 



Since X]i=i ^i converges to a Poisson random variable of parameter s as A ^ 0, it follows that 
'YldLi ^i converges in distribution to a Poisson variable of parameter s, when first 7\/ — > cx3 
and then A ^- 0. On the other hand, using (|4.16p we have that 

^s/A s/A \ 

■ I i=l j 

and therefore, one concludes that also X]i=i ^^ converges in distribution to a Poisson variable 
of parameter s, when first M — > cx) and then A — !► 0. This yields the corollary. D 
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